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Abstract 

Lycurus consists of three species, all with paired, single-flowered spikelets (the 
lower is short pedicellate and usually staminate, the upper long pedicellate and 
perfect). Lycurus setosus occurs in the southwestern USA, northern Mexico, north¬ 
western Argentina, and Bolivia. Allozyme data were used to evaluate genetic di¬ 
versity within and among populations of this amphitropical disjunct species. Elec¬ 
trophoretic examination of 18 putative enzyme loci in 13 populations revealed high 
levels of genetic variation (P ranging from 0.43 to 0.79; H from 0.31 to 0.62) and 
high levels of genetic diversity (F ranging from -0.38 to —1.00). All populations 
possess high levels of heterogeneity ( F is approaching —1, mean of —0.723) and 
exhibit lower levels of genetic fixation among populations ( F ST mean of 0.256). A 
comparison of genetic identity values among populations from North and South 
America indicates that the genetic variation is greater (/ = 0.89) in North America 
than in South America (/ = 0.94), and populations from South America lack six 
alleles found in the North American populations. There was one unique allele found 
in populations from South America. It seems likely that Lycurus setosus has re¬ 
cently dispersed to South America because the populations there contain less ge¬ 
netic variation. 


Resumen 

Lycurus consiste en tres especies, todas con espiguillas dispuestas en pares, flosculo 
unico (el basal brevemente pedicelado y usualmente estaminado, y el distal larga- 
mente pedicelado y perfecta). Lycurus setosus habita en el suroeste de Estatos Unidos, 
norte de Mexico, noreste de Argentina y Bolivia. Mediante el analisis de alozimas 
se evaluo la diversidad genetica dentro y entre poblaciones de esta especie disjunctiva 
anfitropical. El examen electroforetico de 18 loci putativos enzimaticos en 13 poblaci¬ 
ones, revelo altos niveles de variacion genetica (P varfa de 0.43 a 0.79; H varfa de 
0.31 a 0.62) y altos niveles de diversidad genetica (F varfa de —0.38 a —1.00). Todas 
las poblaciones poseen altos niveles de heterogeneidad dentro de las mismas pobla¬ 
ciones ( F IS cerca —1, media de —0.723) y exhiben bajos niveles de fijacion genetica 
entre poblaciones ( F ST media de 0.256). Una comparacion de valores de identidad 
genetica entre poblaciones de Norte y Sur America indican que la variacion genetica 
es mayor (/ = 0.89) en Norte America que en Sur America (/ = 0.94), y poblaciones 
de Sur America carecen de seis alelos encontrados en poblaciones de Norte America. 
Se hallo un unico alelo en poblaciones de Sur Americanas. Probablemente, Lycurus 
setosus ha sido recientemente dipersado hacia Sur America porque esta poblaciones 
contienen menos variacion genetica. 
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Lycurus Kunth consists of three species restricted to the New 
World: L. phalaroides Kunth, L. phleoides Kunth, and L. setosus 
(Nutt.) C. Reeder. The amphitropical disjunct, Lycurus setosus , oc¬ 
curs in the southwestern U.S. and northern Mexico, and again in 
northwestern Argentina and Bolivia. This species was originally de¬ 
scribed by Nuttall (1848) from plants collected in the vicinity of 
Santa Fe, New Mexico, as a distinct genus Pleopogon setosum Nutt. 
Beal (1896) recognized this taxon as Lycurus phleoides Kunth var. 
glaucifolius Beal. Later authors placed this species as a synonym of 
L. phleoides (Hitchcock 1913, 1937, 1939). It was not until C. Reed¬ 
er (1985) revised Lycurus that we came to recognize L. setosus as 
a third species in the genus. 

The genus is characterized by having paired, single-flowered 
spikelets; the lower short, pedicellate and usually staminate, occa¬ 
sionally sterile or perfect, the upper long pedicellate and perfect 
(Peterson et al. 1997). The first glume is 2- or 3-nerved with two 
awns, 3-7 mm long and the lemma is 3-nerved and awned. Lycurus 
setosus can be differentiated from the other two species in the genus 
by having leaf blades terminating in slender seta 0.5-10(12) mm 
long and acute to acuminate ligules 3-10 mm long. 

Upon describing the genus with two species, Kunth (1816) sug¬ 
gested Lycurus was similar to Phleum L. in habit, with affinities to 
Aegopogon Humb. & Bonpl. ex Willd. Based on anatomical char¬ 
acteristics, Sanchez and Rugolo de Agrasar (1986) suggested Ly¬ 
curus be aligned with Aegopogon and Tragus Haller f. (tribe Zoy- 
sieae). Sanchez and Rugolo de Agrasar (1986) also pointed out the 
anatomical similarity between Erioneuron Nash (tribe Eragrosti- 
deae) and Lycurus. Based upon morphological similarities, Mez 
(1921) transferred Muhlenbergia shaffneri E. Foum., considered a 
synonym of Muhlenbergia depauperata Scribn., to Lycurus. Muh¬ 
lenbergia depauperata and M. brevis C. O. Goodd. share many mor¬ 
phological features with Lycurus , most importantly: paired spikelets, 
2-nerved lower glumes with two awns, and 3-nerved lemmas (Pe¬ 
terson and Annable 1991). Pilger (1956) erected the subtribe Ly- 
curinae, which included Lycurus and Pereilema J. Presl. Clayton 
and Renvoize (1986) and Valdes-Reyna and Hatch (1991) took the 
traditional view and place Lycurus in the subtribe Sporobolinae, 
along with Calamovilfa (A. Gray) Hack, ex Scribn. & Southw., 
Crypsis Aiton, Muhlenbergia Schreb., Pereilema , and Sporobolus 
R. Br. Based on chloroplast DNA evidence, Lycurus appears to be 
firmly embedded in the subtribe Muhlenbergiinae, subfamily Chlor- 
idoideae, along with Bealia Scribn., Blepharoneuron Nash, Cha- 
boissaea E. Foum., Muhlenbergia , and Pereilema (Duvall et al. 
1994; Peterson et al. 1995, 1997). 

The origin of the Chloridoideae (=Eragrostoideae, sensu Pilger 
1954, 1956) appears to be in southwestern Africa, where the climate 
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is extremely arid, summer rainfall is common, and mean winter 
temperature is above 10°C (Hartley and Slater 1960). Chloridoid 
centers of diversity occur in southwestern Africa, northcentral (Ti- 
besti, Sahara) Africa, and northwestern USA/northem Mexico 
(Clayton 1975; Hartley and Slater 1960). This wide distributional 
pattern suggests that the subfamily is an old one and that subsequent 
radiation from the African continent has occurred (Hartley and Slat¬ 
er 1960). 

Although there is some doubt as to the specific taxon surveyed, 
chromosome counts for Lycurus setosus indicate that it is tetraploid 
(2 n = 40) with a base chromosome number of x = 10 (Avdulov 
1931; Gould 1964, 1965; J. Reeder 1967, 1971, 1977). The U.S. 
National Herbarium (US) has four specimens of L. setosus from 
Chihuahua and Durango, Mexico, collected by Reeder and Reeder 
(4877, 4884, 4890, 4898) that indicate a chromosome count of 2 n 
= 40. 

The present study, the first analysis of soluble enzymes in Lycu¬ 
rus , was initiated to estimate the genetic diversity within and among 
populations of L. setosus. We also hoped to gain new insights into 
the phytogeographical history and evolutionary processes operating 
among populations with amphitropical disjunct distributions. Similar 
studies of allozyme variation in Bealia mexicana Seribn., Chabois- 
saea atacamensis (Parodi) P. M. Peterson and Annable, C. decum- 
bens (Swallen) Reeder and C. Reeder, C. ligulata. E. Foum., C. 
subbiflora (Hitchc.) Reeder and C. Reeder, Muhlenbergia argentea 
Vasey, M. lucida Swallen, M. torreyi (Kunth) Hitchc. ex Bush, and 
Scleropogon brevifolius Phil, have revealed high intraspecific vari¬ 
ability (H ranging from 0.19 to 1.00) and high levels of genetic 
diversity (F ranging from 0.073 to —1.000) (Peterson and Columbus 
1997; Peterson and Herrera A. 1995; Peterson and Ortiz-Diaz in 
preparation; Peterson et al. 1993). 

Methods 

Thirteen populations representing 365 individuals were sampled 
from throughout the geographic range of Lycurus setosus (Table 1). 
Fresh leaf blades were collected in the field, placed in 3.6 or 5.0 ml 
cryotubes (NUNC), and frozen on site in liquid nitrogen. 

Sample preparation and electrophoresis of enzymes followed the 
general methodology of Morden et al. (1987). Approximately 300 
mg of mature tissue from each plant was homogenized in up to 25 
drops of grinding buffer together with about 50 mg of sea sand to 
enhance disruption of cells. Extracts were absorbed into 2X11 mm 
Whatman filter paper wicks and stored at — 80°C. Electrophoresis 
was conducted in the four gel/electrode buffer systems (L, M, N, T) 
as described in Morden et al. (1987). Starch concentration was mod- 
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Table 1. Field Collections of Lycurus setosus Analyzed by Enzyme Electro¬ 
phoresis. Vouchers deposited at US. 


ARGENTINA. Juyuy: Depto. Humahuaca: 20 km N of Humahuacaon Hwy 9 
towards Tres Cruces, 17 Mar. 1991, Peterson & Annable 11762. Mendoza: Depto. 
Lujan de Cuyo, 12 km SW of Potrerillos, 26 Feb. 1991, Peterson & Annable 11439. 
Salta: Depto. La Poma, 4 km N of Saladillo on Hwy 40, 15 Mar. 1991, Peterson & 
Annable 11746.. Depto. San Carlos, just E of Isonza and 26 km N of Amblayo, 13 
Mar. 1991, Peterson & Annable 11724. San Juan: Depto. Zonda, 37 km SW of 
Zonda at Estancia Maradona, 29 Feb. 1991, Peterson & Annable 11513. Tucuman: 
Depto. Tafi del Valle, 30 km SE of Amaicha del Valle on Hwy 307, 8 Mar. 1991, 
Peterson & Annable 11622. 

MEXICO. Sonora: Sierra El Gato, 5.7 mi E of Huachinera, 9 Oct. 1992, Peterson 
& Annable 12350. 

U.S.A. Arizona: Coconino Co., 6 mi NW of Ash Fork on Canyon Road (39) below 
Antolini, 28 Sep. 1992, Peterson & Annable 12174. Gila Co., 3.7 mi S of Young on 
Hwy 288 towards Globe, 30 Sep. 1992, Peterson & Annable 12218. Pima Co., Santa 
Rita Mountains, 4.8 mi W of Hwy 83 on Forest Service 231 and 0.2 mi W on Forest 
Service 4053, 1 Oct. 1992, Peterson & Annable 12246. Yavapai Co., 8.7 mi SW of 
Jerome on Hwy 89A, 28 Sep. 1992, Peterson & Annable 12184. New Mexico: Grant 
Co., Burro Mountains, 0.3 mi N of Hwy 90 on Mill Canyon Road (Forest Service 
859), 4 Oct. 1992, Peterson & Annable 12299. Hidalgo Co., Peloncillo Mountains, 
10 mi W of Hwy 338 on Forest Service 63 in Clayton Draw, 5 Oct. 1992, Peterson 
& Annable 12320. 


ified to improve gel handling and improve resolution; all gels con¬ 
sisted of Sigma hydrolyzed potato starch at concentrations of 10.6%, 
12.0%, 11.5%, and 12.0% for the L, M, N, and T systems, respec¬ 
tively. For each population, samples from all individuals were in¬ 
cluded together on the same gel. Selected individuals from different 
populations were then analyzed together for purposes of interspecific 
and interpopulational comparisons. 

Gels were sliced and stained for the following 14 enzymes: as¬ 
partate aminotransferase (AAT), aconitase (ACO), adenylate kinase 
(ADK), aminopeptidase (AMP), fructokinase (FRK), glutamate de¬ 
hydrogenase (GDH), glutamate-pyruvate transaminase (GPT), iso¬ 
citrate dehydrogenase (IDH), malate dehydrogenase (MDH), 6- 
phosphogluconate dehydrogenase (PGD), phosphoglucose isomer- 
ase (PGI), phosphoglucomutase (PGM), shikimate dehydrogenase 
(SAD), and triose phosphate isomerase (TPI). Banding pattern in¬ 
terpretations are based on known subunit structure and conserved 
number of loci at the diploid level. Only the faster migrating bands 
of IDH, assumed to be the nuclear-encoded plastid form were sur¬ 
veyed (Gottlieb 1982; Weeden and Wendel 1989). Loci were des¬ 
ignated sequentially with the most anodally-migrating locus desig¬ 
nated 1, the next 2, and so on. Alleles were designated sequentially 
with the most anodally-migrating allele given an a , the next b, and 
so on. 

Values for Nei’s (1972) genetic identity (/) and distance measures 
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Table 2. Genetic Variation in Populations of Lycurus setosus : sample size (n); 
mean number of alleles per locus (A); mean proportion of polymorphic loci (P); 95% 
criterion, mean heterozygosity (H), direct count estimate; and mean fixation index 
(F). 


Collection No. 

n 

A 

P 

H 

F 

11439 

29.0 

1.4 

42.9 

0.429 

- 1.000 

11513 

30.3 

1.9 

71.4 

0.311 

-0.390 

11622 

27.0 

1.8 

57.1 

0.532 

-0.833 

11724 

29.0 

1.8 

57.1 

0.520 

-0.789 

11746 

30.0 

1.7 

64.3 

0.569 

-0.778 

11762 

30.0 

1.8 

64.3 

0.417 

-0.583 

12174 

28.0 

1.9 

57.1 

0.372 

-0.527 

12184 

28.0 

1.8 

57.1 

0.520 

-0.722 

12218 

27.0 

2.1 

78.6 

0.450 

-0.444 

12246 

28.0 

1.9 

78.6 

0.615 

-0.607 

12299 

27.9 

1.7 

57.1 

0.334 

-0.413 

12320 

27.1 

1.9 

78.6 

0.398 

-0.379 

12350 

23.0 

2.0 

71.4 

0.385 

-0.433 

Total Means 

28 

1.8 

64 

0.45 

-0.61 


were computed for pairwise comparisons using BIOSYS-1 (Swof- 
ford and Selander 1989). Although all these populations are prob¬ 
ably tetraploid in origin, an initial BIOSYS-1 run was undertaken 
to measure each population’s genetic variability and to compute ge¬ 
netic identities (Nei 1972). Standard measures of genetic variation 
(Table 2) were computed including mean number of alleles per locus 
(A), proportion of polymorphic loci (P), mean heterozygosity ( H ), 
and mean fixation index ( E) which measures the deviation of ge¬ 
notypic proportions from Hardy-Weinberg expectations (Wright 
1965). The distribution of genetic variation was determined using 
F-statistics where F IS is the fixation index within populations, F IT is 
the overall fixation index or inbreeding coefficient, and F ST measures 
the degree of differentiation among populations (Wright 1965, 1969; 
Jain and Workman 1967). The patristic distance matrix was calcu¬ 
lated using the Prevosti distance index (Wright 1978) and after op¬ 
timization of branch lengths, a corresponding phenogram (Fig. 1) 
using the Wagner procedure was produced (Swofford and Selander 
1989). 


Results 

Eleven enzyme systems encoded by 18 putative loci were con¬ 
sistently scorable by starch gel electrophoresis: AAT-1, AAT-2, 
ACO-1, ACO-2, AMP-1, AMP-2, GDH, GPT, IDH, PGD, PGI-1, 
PGI-2, PGM-1, PGM-2, SAD-1, SAD-2, TPI-1, TPI-2. Several en¬ 
zymes or putative loci, viz., ADK, FRK, MDH, and PGD-2 were 
not scorable due to faint or inconsistent staining. Allele frequencies 
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Fig. 1 . Phenogram showing genetic distance among populations of Lycurus setosus. 
Correlation coefficient = 0.841; length = 0.971; numbers refer to population collec¬ 
tions given in Table 1; populations from South America are marked with an asterisk; 
scale indicates distance from midpoint. 


for all 13 populations surveyed appear in Appendix 1. The number 
of alleles per polymorphic locus ranged from two in AAT-1, AAT- 
2, AMP-1, AMP-2, GDH, GPT, PGD, PGM-1, PGM-2, TPI-1, and 
TPI-2 to six in PGI-1. The following loci were fixed for a single 
allele for all populations: AMP-1 and AMP-2. The largest number 
of alleles per locus in a population was four, occurring in PGI-1 
(population 11513). There was no evidence of duplicated loci, how¬ 
ever, AAT-1, AAT-2, GDH, and GPT were fixed for a pair of dif¬ 
ferent alleles and removed from the genetic analysis. 

The sample size (n) per population ranged from 23 to 30 individ¬ 
uals and the mean number of alleles per locus (A) within populations 
ranged from 1.4 to 2.1 (Table 2). The percentage of polymorphic 
loci ( P ) ranged from 0.429 to 0.786 and the mean heterozygosity 
(//), direct count estimate, ranged from 0.311 to 0.615, indicating a 
high level of heterozygosity at most polymorphic loci (Table 2). The 
mean fixation index (F) within populations, or inbreeding coeffi¬ 
cient, ranged from —0.38 to —1.00. Genetic variability as measured 
in pooled populations from North and South America was: A = 1.7, 
1.9; P — 60, 68; H = 0.46, 0.44 (Table 3). 

There were 45 alleles recorded in the North American populations 
and 40 in the South American populations (Appendix 1, Table 3). 
One unique allele was detected in each of four populations 11513 
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Table 3. Genetic Variability in North and South American Populations of 
Lycurus setosus : number of populations sampled (n); mean number of alleles per 
locus (A); mean proportion of polymorphic loci (P); 95% criterion, mean heterozy¬ 
gosity (H), direct count estimate; and alleles in common and unique (u). 


Country 

n 

A 

P 

H 

Alleles (u) 

South America 

6 

1.7 

59.5 

0.463 

40 (1) 

North America 

7 

1.9 

68.4 

0.439 

45 (6) 


(PGI-l-b), 12184 (PGI-2-c), 12320 (TPI-l-a), and 12350 (SAD-1- 
a). Two of these unique alleles were found in single individuals in 
populations 11513 and 12320; PGI-2-c (12184) was found in two 
individuals and SAD-l-a (12350) was found in eight individuals. In 
North America all seven populations shared alleles PGI-l-a & d and 
two populations (12174, 12350) shared allele SAD-2-a. None of 
these three alleles was present in South American populations. In 
summary, a total of six alleles are exclusively shared by North 
American populations and a one allele is unique in South America 
(Table 3). 

Partitioning of genetic diversity, or the fixation of alleles at dif¬ 
ferent hierarchical levels, within and among populations of Lycurus 
setosus was calculated using F-statistics where the fixation index 
within populations (F IS ) ranged from —0.019 to —0.927 (mean 
—0.723). The amount of genetic diversity among populations (F ST ) 
ranged from 0.017 to 0.456 (mean 0.256). The overall fixation index 
(F it ) ranged from 0.008 to —0.758 (mean —0.281). The primary 
component of F IT was F /5 , i.e., the F ST values were much smaller, 
indicating greater heterogeneity within populations than among 
them. 

A phenogram (Fig. 1) summarizes the interpopulational relation¬ 
ships based on genetic distance values. Six populations from South 
America are differentiated on a single branch at the 0.38 level or 
distance from root, whereas seven populations from North America 
are differentiated on a single branch at the 0.20 level. These two 
branches separating the North from the South American populations 
are then joined near the base at the 0.009 level. 

Mean genetic identities (Nei 1972) among populations were quite 
variable ranging from 0.718 between population 12184 and 12350 
(both from North America), to 0.978 between population 12299 and 
12320 (both from North America). The mean genetic identity for 
all 13 populations was 0.89. A comparison of identity values among 
populations from North and South America indicates that the genetic 
variation is greater (/ = 0.89) in North America than in South Amer¬ 
ica (/ = 0.94). 
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Discussion 

All populations of Lycurus setosus examined in this study show 
high levels of genetic variation comparable to that found in out- 
crossing plant species (Hamrick and Godt 1989; Hamrick et al. 
1979). This is reflected in the high mean values for all populations 
for proportion of polymorphic loci of 0.64, heterozygosity of 0.45, 
and negative F values of —0.61. All of these values indicate a con¬ 
sistent excess of heterozygotes (Table 2). All populations possess 
high levels of heterogeneity within populations (F IS approaching — 1, 
mean of —0.723) and exhibit lower levels of genetic fixation among 
populations ( F ST mean of 0.256). 

Since the phenogram (Fig. 1) was computed using genetic dis¬ 
tance it mirrors the results inferred from genetic identities. Clearly, 
populations from South America form a more compact group re¬ 
flecting a lower genetic variability. However, population 11513 from 
San Juan, Argentina, is quite variable (mean number of alleles per 
locus, A = 1.9, and proportion of polymorphic loci, P = 71.4) and 
similar to North American populations. 

Genetic identities values indicated that populations from South 
America (/ = 0.94) were more similar to each other than those from 
North America (/ = 0.89). Both of these values fall within the range 
of other intraspecific identity values obtained for plant populations 
and are considerably higher than those reported for congeneric spe¬ 
cies (Gottlieb 1981 ,1 = 0.67). Other intraspecific identity values (/) 
for eragrostoid grasses reported are: 0.96 for Bealia mexicana ; 0.94 
for Chaboissaea atacamensis\ 0.82 for C. ligulata ; 0.88 for C. sub¬ 
biflora ; 0.95 for Muhlenbergia argentea ; 0.98 for M. lucida\ and 
0.93 for Scleropogon brevifolius (Peterson and Columbus 1997; Pe¬ 
terson and Herrera A. 1995; Peterson et al. 1993). Because there 
was a single unique allele present in the South American populations 
where less genetic variation exits and these same populations were 
missing six alleles found in North American populations, it seems 
likely that Lycurus setosus has recently dispersed to the Southern 
hemisphere. 

Scleropogon brevifolius (Peterson and Columbus 1997) and Muh¬ 
lenbergia torreyi (Kunth) Hitchc. (Peterson and Ortiz, in preparation) 
apparently have similar biogeographical histories. Like Lycurus seto¬ 
sus , populations of both species have lower levels of genetic variation 
in the Southern hemisphere and have very few or no unique alleles 
while lacking alleles shared by North American populations. This 
North to South American dispersal is seen in Chaboissaea , where a 
single species or vicariad, C. atacamensis , has recently dispersed from 
Mexico to South America from its closest sister, C. ligulata (Sykes et 
al. 1997). The north to south migration pattern is not always the case, 
and the reverse (south to north) is exhibited in Bothriochloa Kuntze 
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(Allred 1981) and Erioneuron Nash (Peterson in preparation). Cha- 
boissaea, Lycurus, and Muhlenbergici, all members of subtribe Muhl- 
enbergiinae, have centers of diversity in North America. In contrast, 
Erioneuron , a member of the Munroinae (including Blepharidachne 
Hack., Dasyochloa Willd. ex Rydb., and Munroa Torn) has a center 
of diversity in South America. Two major themes that all of these 
disjunctions seem to have in common is 1) recent migration, since little 
morphological differences and very little genetic differentiation exist, 
and 2) the putative dispersed species or taxon has migrated away from 
a center of diversity. 
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Appendix 1. Allele frequencies for 13 populations of Lycurus setosus (see Table 1 for population localities) from North and South America. 
Unique alleles found in a single individual are marked with an asterisk. 
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Appendix 1. Continued. 
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